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Abstract. Anderson localization is known to be inevitable in one dimension for 
generic disordered models. Since localization leads to Poissonian energy level statistics, 
we ask if localized systems possess “additional” integrals of motion as well, so as 
to enhance the analogy with quantum integrable systems. We answer this in the 
affirmative in the present work. We construct a set of nontrivial integrals of motion for 
Anderson localized models, in terms of the original creation and annihilation operators. 
These are found as a power series in the hopping parameter. The recently found 
Type-1 Hamiltonians, which are known to be quantum integrable in a precise sense, 
motivate our construction. We note that these models can be viewed as disordered 
electron models with infinite-range hopping, where a similar series truncates at the 
linear order. We show that despite the infinite range hopping, all states but one are 
localized. We also study the conservation laws for the disorder free Aubry-Andre 
model, where the states are either localized or extended, depending on the strength 
of a coupling constant. We formulate a specific procedure for averaging over disorder, 
in order to examine the convergence of the power series. Using this procedure in the 
Aubry-Andre model, we show that integrals of motion given by our construction are 
well-defined in localized phase, but not so in the extended phase. Finally, we also 
obtain the integrals of motion for a model with interactions to lowest order in the 
interaction. 
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1. Introduction 

The simplest theoretical model to study localization for non-interacting particles in the 
presence of disorder was proposed by Anderson [T] . A single particle localized state has 
a wavefunction that decays exponentially abont some point in space over a characteristic 
localization length. In three dimensions, localized states exist below a certain energy 
(the mobility edge) for a given strength of disorder. A disordered electronic system is 
thus localized if its Fermi energy lies below the mobility edge. In one and two dimensions, 
an infinitesimal amount of disorder is sufficient to localize all single particle states and 
thus a disordered non-interacting electronic system is always localized [2113!. 

Recent developments in the area of eigenstate thermalization mu relate closely 
to the above well established notions of Anderson localization. In this context, it is 
believed that an isolated localized eigenstate does not thermalize, in the sense that no 
subsystem of it can be brought into thermal equilibrium by exchanging heat with the 
rest of the system. An analogons statement can be made about information, as defined 
through an appropriate partial trace of the density matrix. A related feature of such a 
system is the lack of level repulsion in its energy level spectrum. This can be thought 
of as arising from the presence of almost degenerate states localized so far apart that 
they are unable to hybridize to lift the degeneracy. 

The effect of interactions on snch systems is very interesting. Interactions among 
the elementary degrees of freedom generically tend to drive the system towards 
thermalization and delocalization laiH]. This tendency competes with the the one that 
causes localization in the presence of disorder. Understanding the resultant phenomenon 
of many body localization, that is observed for sufficiently strong disorder, is currently 
a very active area of research punmiiiniiisi. 

Another class of systems that fail to thermalize are integrable ones. These often 
contain a variable parameter (snch as interaction or external field strength, which 
we denote here as y) and possess a set of similarly dynamical (i.e. depending on 
the parameter) integrals of motion. Standard examples of snch systems are the one¬ 
dimensional Hubbard and XXZ models. In these examples, the integrals of motion A 
are polynomial in y with the order of the polynomial [HI (HI [lEl HZl UHl UHl ES El] 
increasing with k. An arbitrary linear snperposition of all integrals Q = A “ also 

an integral in its own right - is an infinite power series in y. Gaudin magnets [2211^ on 
the other hand provide examples of integrable models where all conserved chargeqil are 
linear vn y. It should be emphasized that there is no generally accepted precise notion 
of integrability in quantum mechanics [2H [25] in contrast to classical mechanics where 
it is unambigous. However, we do not dwell on this issne in present worli^. The only 
aspect important for us here is the existence of parameter-dependent conservation laws. 

t We use the term conserved charges interchangeably with conservation laws or integrals of motion. 

§ By integrable we will generally mean quantum many-body models colloquially recognized as such, 
see examples in this paragraph. The only exception are type-1 Hamiltonians that stem from a recently 
proposed well-defined notion of quantum integrability |26l |2Il EH ESI ESI [30] . 
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Conserved charges greatly constrain the dynamics of integrable systems. As a 
result, when started off from an arbitrary initial state in isolation, these systems do not 
evolve in a way that causes thermalization in the sense of the above paragraph urn- 
Additionally the usual space time symmetries result in degeneracies in the energy level 
spectrum, and hence a lack of level repulsion [32]. The addition of perturbations 
destroys such conservation laws and restores level repulsion, although the strength of 
the perturbations has a non-trivial finite-size dependence [3Sl EH EH] • 

In this context, it is natural to ask in what ways are localized systems similar to 
integrable ones. In particular we may ask if (parameter dependent) conservation laws, 
similar to those in integrable systems exist for localized systems. It has been argued in 
the context of many-body localization that they do, and results related to the growth of 
entanglement in these systems are predicated on their existence [36] IHT] [38]. However, 
obtaining the structure of the conserved charges directly in terms of microscopic 
parameters remains a challenge and effective renormalization procedures need to be 
employed instead I39i unj. The situation is less complicated in the absence of 
interactions since the Hamiltonian is that of a single particle. Nevertheless, obtaining 
the charges systematically and analytically in terms of the microscopic parameters of 
the Hamiltonian is non-trivial. In this paper we outline the procedure to do so. We also 
elucidate the connection between localization and conserved charges. 

In this work we study a general one dimensional model with on-site disorder that can 
interpolate between models with long-range hopping and the more standard Anderson- 
type one. The starting point is a Type-1 Hamiltonian reviewed in Refs. [SHlETlEHl 
EH ESI ED). This was introduced as the most simple model of quantum integrability 
in finite dimensional spaces. This model has infinite ranged hopping, and as such has 
no inbuilt metric or length scale. We first show by calculating its Participation ratio 
(PR) [ini E] the perhaps surprising result that all states except one are localized. This 
is done as follows: an eigenstate If)) of the Hamiltonian is expanded in a basis of position 
eigenstates on the lattice as |'0) = ^k\k), where k labels the position eigenstates and 

Ck are the coefficients in the expansion. The PR for this state is then defined as 


PR^ = 


iEk\ck\^f 

Ek\ck\^ 


( 1 ) 


It is usually understood that PR^ ^ 0(1) indicates localization while PR^ r\j 0(N)- 
delocalization, where N is the number of sites. While this definition is valid for a fixed 
wave function, we may also define the PR at a given energy, as later in the paper, where 
an averaging over disorder realizations is carried out, at a fixed energy. 

The Type-1 model has a known set of conservation laws, which inspire the 
construction of a generic Anderson-type model having only nearest-neighbor hopping. 
In Id it is well known that for this model, all single particle eigenstates are localized for 
any strength of the disorder. The conserved charges of this model are then constructed 
by analogy with the Type-1 Hamiltonian. These charges are expressed as a power series 
in the hopping, whose coefficients we determine by means of an algorithm. We also 
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show that the series, upon disorder averaging over a “non-resonant” ensemble- defined 
below, is convergent. This provides numerical evidence that the ensemble chosen and 
the procedure of averaging the coefficients in the conserved charges over the ensemble 
is meaningful. 

We then turn our focus to a model which contains both localized and delocalized 
phases (i.e. phases in which a//single particles states are either localized or delocalized). 
This is the Aubry-Andre model in which the random potential is replaced by a 
quasi-periodic one. This allows us to test our criterion for the convergence of the power 
series and clearly elucidate the connection of the conserved charges to localization. Thus, 
the convergence (divergence) of the power series representation of conserved charges can 
indeed be identified with the presence (or absence) of localization and the localization- 
delocalization transition can be located using the charges. Finally, we investigate the 
effect of interactions and argue that a power series in the interaction becomes intractable 
and thus obtain the the conserved charges only to first order in it. 

We emphasize that the main feature of our construction is that the conservation laws 
do not depend explicitly on the wavefunctions of the single particle energy eigenstates. 
In fact, the recursion relations we obtain for the coefficients of the expansions of the 
conserved charges are the same for all generic one dimensional models. Our approach 
is thus completely model independent requiring no knowledge of exact solutions or 
properties of energy eigenfunctions. 

Another important aspect of the construction of conservation laws we emphasize 
here, which has not been discussed before is ‘gauge freedom’ of a certain kind, defined 
more precisely later. We show that a judicious choice of gauge can bring out important 
features of the conserved charges, such as the truncation of their series representation at 
finite order. These features can be obscured in gauges that arise in constructions of these 
charges from direct applications of standard methods such as the Rayleigh-Schrodinger 
series or the locator expansion. 

2. Lattice Models 

We consider a general Hamiltonian of non-interacting particles hopping on a one 
dimensional lattice with an on-site potential 

H = H{y) = Y^ CiHi Ujc\cj, (2) 

i ij 

where ct and q are fermionic creation and annihilation operators, n, = ctcj is the number 
operator, Cj is the on-site disordered potential, and Uj is the hopping between sites i 
and j. The parameter y is a, real number introduced for convenience, which; it allows 
us to perform an expansion of the conserved charges in its powers. 

Our general strategy to construct construct conserved charges for this models will 
be to first consider the ‘unperturbed’ Hamiltonian which only has the on-site potential. 
The conserved charges for this Hamiltonian are simply the operators n*, which are 
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independent and commute with each other and the Hamiltonian. It can also be readily 
seen that the eigenstates of this Hamiltonian are completely localized on the individual 
sites. Thus the zeroth order Hamiltonian trivially describes a localized system with 
conserved charges. We now show that upon introducing the hopping, new conserved 
charges Qi appear, which can still be labeled by the site indices i while the system 
remains localized. To do this, we consider different types of hopping parameters tij. 


3. Type-1 Hamiltonians 


We now summarize a known set of conserved charges Qj. We rework the construction 
in Refs. [SZIEHIES], in a fashion that suggests a natural generalization for short ranged 
models. These charges are linear in the hopping (or the parameter y), and commute 
exactly with the Hamiltonian of the Type-1 family. The Type-1 Hamiltonian is obtained 
from Eq. (|2]) by specializing to infinite ranged hopping tij = with arbitrary 

parameters jj. Specializing to j = 0 we write down the charge Qq 


Qo = no-y 






^0 ~ 


[tokiclck + c\co) -alno - 


(3) 


where q;° and /5° are yet to be determined.The commutator of Qq and H vanishes to 
linear order in y by construction.The surviving term is of 0{y‘^) and is given by 

[<3o, —-— [A{0,j, k){clcj - c]co) 

jk 

+ k){clcj - c]ck)] = 0, (4) 

where 


^Ok^kj ^k^Oj 
k) = tojtok - Pktkj- 

(5) 

A few words on the form of Eq. (El) are appropriate here. The last two terms 
—a^Hk —/dfcUo commute with H{y = 0) trivially, since they are expressed in terms of the 
number operators. These actually represent a particularly convenient “gauge choice”, 
their presence enables the second order term 0{y‘^) to vanish, and thus the commutator 
series to truncate exactly for the Type-1 matrices. The requirement that [Qo, W] = 0 is 
satisfied by the following form of tij. 


tij — Ailj 



( 6 ) 


this gives A{f),j,k) = B(0,j,k) = 0. It is straightforward to extend this definition to 
arbitrary Qj, and further to show that [Qi, Qj] = 0 Vi, j, so the operators Qi are indeed 






Integrals of motion for one-dimensional Anderson localized systems 


6 


the conserved charges of the Hamiltonian i/[2Sl EZl- The Hamiltonians described by 
tij of the form given in Eq. (jHl) are called Type 1 (23 ES], and can also be interpreted 
geometrically as representing a ‘d-simplex ’ H- 


3.1. PR for Type-1 Hamiltonians 


All single particle states of Type-1 Hamiltonians (El), except possibly the ground state 
for y > 0 or the highest energy state for y < 0 are localized, see e.g. Fig. [3 

This can be understood in more detail from the exact solution for the spectrum 
of these models j2Z]- Exact un-normalized single particle eigenstates of the Hamiltonian 
(Ej) are 

N-l f 
i=0 ^ 

and the corresponding eigenvalues E (energies) are solutions of the equation 


N-l 


E 


t? 


E - ti 



( 8 ) 


Suppose Cj are ordered in the ascending order. By plotting the left hand side of Eq. (jH]) 
as a function of E, one can verify that it has N — 1 real roots Ei,E 2 ,... E^v-i located 
between consecutive Cj, i.e. Cj-i < Ei < ei. The remaining root Eq is also real and is 
below Co (ground state) for y > 0 and above ev-i for r/ < 0 (highest excited state). 

Eqs. d?]) and dH]) also provide an exact solution for one fermion (Cooper) pair and 
one spin flip sectors of the BCS and Gaudin models, respectively. 


Hbcs ^ y ^ 

Hi{x) = si -yY^ , 

. / . Q 


(9) 


where are spin-full fermions and Sj are quantum spins of arbitrary magnitudes s*, see 

Ref. [27] for details. For the BCS [Gaudin] model one needs to replace c\ 

Eq. dZ])) 7i 1 [v^2ii] and the corresponding eigenvalue is equal to 2E [2sj(£' —ej)“^] 
rather than E. Our results for the PR of Type-1 Hamiltonians therefore also apply to 
these sectors of these models. 

The PR defined through Eq. (ITTl) reads 


PR£; = 




(E-eX 


E 



(E - £.)* 


( 10 ) 


For concreteness we take y > 0. Then, the ground state is Eq < cq. We assume that most 
7 j are of the same order of magnitude and consequently the vector with components y* 
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is delocalized. Further, we take e* to lie in a fixed interval that does not scale with N, 
e.g. from —w to w. 

For excited states Ek is between Ck-i and e^. The summations in the numerator 
and denominator of Eq. TO both come from Cj in a small vicinity of for large N and 
converge as respectively, where n = \i — k\. The numerator and 

the denominator scale as and 7^/h^, where 5 oc l/A/” is the mean level spacing 

between e* in the vicinity of e^. Therefore, PRs^, is of order 1 (much smaller than N) 
meaning excited states are always localized. Fig. H] shows PR for N = 10^ uncorrelated 
random e* uniformly drawn from an interval (—1,1) and the same distribution of q*. 

Consistent with our numerical results, we estimate the largest PR for excited states 
to scale as InN, i.e. 

PR^ff^alnA^, (11) 

for large N, where a depends on N much weaker than In N. Such values of PR come 
from clustering in Cj. Indeed, suppose spacings h* = Cj+i — e* between m of Cj for i 
from k to k m are all much smaller than 5k-i and, moreover, Ck+m — < ^k-i- R 
follows from Eq. TO that PRej, ~ PREfc+„^+i ~ ni because the above e* contribute 
most to these PRs. Normalized spacings Sj = 5i/5 are distributed according to the 
Poisson distribution P{s)ds = e~^ds. The probability of having m spacings between 0 
and So 1 is then roughly s™. We need msQ < 1 and also Ns]f = 1 so that at least 
one such clustering occur4||l . This implies m ~ In A^/ln(ln A^) and Eq. flTTl) follows. 
Numerically we find that typical values of a ~ 1 — 3 and averaged over disorder a ~ 1.7, 
at least for N” = 2^ — 2^^. Note that according to this argument such large values of PR 
typically come in pairs spaced by m + 1, roughly equal to the value of the PR itself. We 
also stress that, in contrast to the largest PR, a typical (and average) PR is something 
between one and three for any N (does not scale) as can be seen from Pig. [H 

It is interesting to compare this In N behavior to the fiat band localization studied 
earlier [M Ii5] . The latter leads to a (weakly) divergent PR in the localized regime, a 
phenomenon that is viewed as corresponding to critical (power law type) localization. 
The Type-1 Hamiltonian kinetic energy may also be viewed as a “flat band” model, with 
a flat dispersion for all except one state. Indeed, for tij = ji'-jj all but one eigenvalues of 
the second term in Eq. (E]) are zero. The non-zero eigenvalue (ground state for y >0) 
corresponds to the eigenstate qjc||0). 

Let us consider limits y ^ 0 and y ^ oo separately. When y ^ 0 all states 
are localized as expected. Indeed, Eq. dHj) implies Ek e^, summations in Eq. fITOl) 
are dominated by the i = k term and we obtain PRe^. = 1 for all k. When y ^ oo 
excited states are localized as before because Ek for /c > 1 remains trapped in the 
interval (e^-i, e^). The ground state energy on the other hand diverges - Eq. (|H)) implies 

II More precisely, the probability that m of occur in an interval of length S for Poisson distribution 
is which however still leads to the same estimate m 
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Eq — 7 - Then, €j are negligible as compared to Eq in Eq. (flOl) and 


PReo 


T.iit 


( 12 ) 


which is of order N according to our choice of 7 j. The ground state is therefore 
delocalized for y —>■ oo. It undergoes a localization-delocalization crossover at a certain 
yc, which we estimate below in this section. 

It is possible to evaluate the PR analytically to leading order in 1 /A/^ for distributions 
of Cj and 7 * with negligible short range fluctuations (such that the spacing h* = Cj+i — e* 
changes slowly with z - is of order 1/N for all i - and similarly for 7 j). For 

simplicity, let us take constant 7 j, which we can set to one with no loss of generality, 
and equally spaced e*, i.e. 5i = 5 = 2w/N. 

For excited states, we write E^ = Ck — akS, where 0 < < 1, and solve Eq. (|H1) for 

to the leading order in 1 /A^ as described in Appendix B of[lEj. This yields 


cotTrafc = —-|--ln ^^~*~^ = /(cfc)- (13) 

TVy TV W - Ck 

We note that X = y/S is the proper dimensionless coupling constant in the sense that it 
must stay flnite in the N ^ 00 limit. This is because the second summation in Eq. ((21) 
scales as N‘^ for tij = 7 j 7 j and our choice of 7 ^. Therefore, we need r/ochocl/A^so 
that both terms in Eq. ((21) are extensive in the thermodynamic limit. For the BCS 
Hamiltonian in Eq. ((91), so defined A is the dimensionless superconducting coupling (TTj. 
Eq. ((Til becomes to leading order in 1/N 




E oo 
n=0 


n=0 I (n+afc)2 


+ 


1 2 


E oo 
n=0 


{n+ak)^ {n+l-akY 


(14) 


which evaluates to 


PRsfc 


3 

1-1-2 COS^ TTftfc 


3 + 3/^(efc) 

1 -|- 3/^(efc) 


(15) 


This answer is in good agreement with numerics already for N = 20, see also Fig. [K 
Note that 1 < PRe^ < 3. 

We saw above that the ground state energy Eq ^ —00 ns y ^ 00 , while Eq —?• cq 
for y ^ 0. Let y be large enough that Eq is well separated from cq. Then, we can 
replace summation in Eq. ((HI) with integration and obtain 


Eq-w S 

in —-= — 

Eq-Iw y 


2w 

W 


(16) 


Performing the same replacement in Eq. f(T(l and using Eq. f(T 6 l) . we derive 


PReo 


3N 

1 - 1-2 cosh(5/|/) ’ 


(17) 
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Figure 1. PR of eigenstates of a Type-1 Hamiltonian for y = .004, N = 10^ in 
ascending order according to the energy. Each and is an independent random 
number uniformly distributed in an interval (—1,1). Larger circle near the left top 
corner indicates the ground state, which is extended. Left inset is the same as above, 
but averaged over 10^ realizations of disorder and compared to Eq. (USD for the same 
w. The right inset shows the PR (except the ground state) for N = 10^ equally 
spaced e^, = 1 and y = .004 similarly compared to Eq. ([T5]) (the two curves are 

indistinguishable). 
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Note that in the limit y —>■ oo, PR£;o = V in agreement with Eq. (fT2l) . This expression 
also allows us to estimate the value Pc beyond which the ground state becomes extended. 
We obtain Ac = Pc/^ ~ 1/ In V. This also corresponds to the coupling for which the 
gap in the spectrum A = Ei — Eq ^ —w — Eq becomes comparable to the spacing 5. 
For a superconductor described by the BCS model (ED this localized-extended crossover 
translates into a normal-superconducting one j4^l49] . As N ^ oo this crossover becomes 
a quantum phase transition at A = 0, i.e. any infinitesimal coupling is sufficient to make 
the ground state extended (superconducting). The localized character of the excited 
states for the specific case of 7 j = 1 has been demonstrated in a previous work as 
well ji3] . 

4. A model with finite-ranged hopping 

We now consider the following Anderson-type model in one dimension with nearest 
neighbor hopping. 

H = eiHi - J^(c|cj+i -F h.c.) 

i i 

= Ho + yH,. (18) 

This corresponds to the case with y = t for \i — j\ = 1 and 0 otherwise for the general 
Hamiltonian in Eq. (ED. Hq is the zeroth order Hamiltonian with only the on-site 
potential and Hi contains the hopping. It is known that all single particle eigenstates 
of this Hamiltonian are localized BEI. 
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Proceeding as for the case of Type-1 Hamiltonians, we focus on the conserved charge 
Qq, corresponding to the site i = 0, which to lowest order is equal to no- However, in 
this case Qq is not simply linear in y. In fact, it can be argued that the an expansion 
of Qo in the hopping does not truncate at any finite order in the thermodynamic limit. 
Indeed, as explained in the Introduction, conserved charges are generally infinite power 
series in y. We thus assume Qi of the form 

Qi = PiO + yPil + y^Pi 2 + • • •, (19) 

where Pjo = and Pn , Pj 2 • • • are operators to be determined in terms of the 
microscopic parameters subject to the condition [Qi,H] = 0. For concreteness, we 
first take our one dimensional system to be a finite-sized ring of Ai -|- 1 sites going from 
0 to N. 

Since the Hamiltonian H and and all the zero order charges n* are quadratic in the 
creation and annihilation operators, we take all the operators Pa, Pi 2 ■ ■ ■ to be similarly 
quadratic, i.e. 

Pim = ^Vj^\i)clcj, (20) 

jk 

where the symmetric coefficients are to be determined. We have 

[Qa H] = [Po, Po] + 5] {[Pim, Pi] + [P™+i, Poj). 

m 

The requirement that the commutator vanishes to all orders in y requires 

[Pim, Pi] + [Pim+l, Pq] = 0, (21) 

and yields a recursion relation among y's. 

( 22 ) 

^a- ^ ^ 

J 

with initial conditions ri^),\i) = diadib- The diagonal term represents a “gauge” 

freedom, since the corresponding term in Pj^ commutes trivially with Pq. We further 
discuss this freedom below. Specializing to the case of nearest neighbor hopping Eq. ca 
and with i = 0, it can be verified that terms present in Pq^ are of the form 

• (CoCm + dco) 

• ^0,lV-(m-l)['^oCAr-(m-l) + cjy_(^_^)Co] 

• Eqym,K-i|=even<m4"'^(4ci + c]q) (if ^ is eveu) 

• Eqym,|i-i|=odd<m4"'^(cl9 + c]ci) (if 771 is odd) 
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Figure 2. Schematic diagram showing hopping terms present in the operators P 1 —P 4 . 
The base site 0 is in the middle and its neighbors are sites 1 and N, since we imposed 
periodic boundary conditions. Lines connecting pairs of sites indicate the presence 
of the corresponding hopping term in the operator Pm- Note that the range of the 
hopping in Pm increases with m. 


This is shown schematically in Fig. [2] for the first few Pm- 

The Qfs are related to each other by translating all site indices in the above relations 
by an appropriate number. By construction, they all commute with H. Since H is 
generally non-degenerate, this implies Qj also commute among themselves, [Qi, Qf = 0 
V z,j. To see this, first recall that for Hermitian matrices [^, B] = [^,(7] = 0 implies 
[5, C] = 0 as long as eigenvalues of A are non-degenerate. All operators involved in 
the above construction of Qi are of the form A = "Yhij ) where Aij is a Hermitian 

N X N matrix, which represents operator A in the sector with total particle number 
n = 1. Moreover, the commutativity of any two such operators is equivalent to that 
of the underlying matrices. Eigenvalues of the Hamiltonian in the n = 1 sector at 
y = 0 are e*, which are assumed to be distinct, i.e. the corresponding matrix is non¬ 
degenerate at y = 0. By continuity of the eigenvalues in r/, it remains non-degenerate 
in some finite interval (until the first level crossing) of the real axis containing y = 0. 
Thus, [Qi,Qj] = Oy i,j in this interval of y. But, as can be seen e.g. from the above 
construction of Qq, commutativity of Qi on any finite interval of values of y implies that 
they commute for all y. 

We noted above that commutation relations (l2Tl) and consequently recursion 
relations fl22l) do not constrain the diagonal part of the coefficients i.e. rIp\ 

for m > 1. The choice of R^\i) however does affect the off-diagonal part of for 
k > m. In our construction of Qi we set R^f^\i) = 0 for all m > 1, since this leads 
to the most compact description of these objects, we will refer to this as the standard 
gauge. Conserved charges Qi resulting from any other choice R^^ {i) uniquely relate to 
our standard gauge Qfs, a brief calculation shows their relationship is 

Qi = Q^ + (23) 

m r 

Another advantage of our choice of a gauge is in a simple relationship between the 
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H = (24) 

i 

To see this, consider the difference 

H-Y. + ..., (25) 

i 

where Wi are ^-independent operators. Note that the zeroth order term cancels in the 
difference. Since H commutes with all Qj, the right hand side (RHS) of Eq. (12^ must 
also commute. This implies in particular [Wi , Uj] = 0 for all i (from the coefficient at 
the lowest power of y in the commutator of the RHS with Qi), which in turn means 
that Wi = Now note that the left hand side (LHS) has zero diagonal matrix 

elements, i.e. no terms of the form cjc^. This is becanse the zeroth order term is absent, 
while higher order terms have no diagonal matrix elements since = 0 for all m > 1 
in onr gange (and similarly the diagonal is absent in other Qi). Then, the diagonal 
matrix elements must vanish on the RHS as well, to all orders in y. In particular, 
r] = 0, i.e. H^i = 0 and 

H-Y = y^^^ + + • • • (26) 

i 

Applying the same argument to the RHS of this equation we similarly obtain W 2 = 0 
etc., until we finally arrive at Eq. (125). 

4 . 2 . Type-1 Hamiltonians redux 

We have seen above that the conserved charges are power series in the hopping. This 
is nnlike the case of Type-1 Hamiltonians, where the power series truncates after the 
first term. The gauge where the series truncates corresponds to having distinct terms 
for m = 1, one can see in Eq. (|3]) (the gauge terms are indicated in the lower braces). 

It is an amusing exercise to determine the correct gauge terms that lead 

to truncation, starting from the recursion relations Eq. (1221) . To obtain Type-1 
Hamiltonians we set Uj — 7i7j, so that the recursions simplify to 

(qrb) - 

br’w (27) 

3 

With the initial condition = diadu,, we obtain at the first level 

oil(i) = (*) + ~ (At - t) ■ (28) 

~ Cfe 

At this point we pause and ask if we can choose the gauge term R^\i) such that ri^i^(i) 

can be made to vanish identically, so that the iterations stop at the first level. Erom 
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Eq. fl?7l) we see that the relevant condition is the vanishing of (Y^^\i) — j • Using 

Eq. (EZD compute 


U'l') = + -f-a - 1) 



(29) 


We may choose so that the term in braces vanishes, thus leading to the truncation 
of the iterations. From Eq. (I28l) we have the complete first order term, and we can 
proceed to construct the charge (denoting the currents by the symbol Q) 

Qi = ni + yY^ 7]^^ {i)clcb, (30) 

ab 


which is identical to that in Eq. (jSI). 

The use of the gauge term here is very special, and guided by our understanding 
of this model. On the other hand, we could by default set all the gauge terms to 
zero, giving us the irreducible (i.e. standard gauge) currents. These no longer truncate 
even for Type-1 Hamiltonians. For completeness we note the second order term for the 
current in this (standard) gauge 


Qi = Ui {y + y^Y^ —^) X (4o + cJq) 

6j — €i €j — 6i 

3 


(31) 


Thus the Type-1 Hamiltonians allow for variety of expressions of the constants of motion. 
To establish their equivalence in general is a subtle problem, where some surprising 
results have been found quite recently in Ref. EDI. 

This type of gauge choice, made explicit in our construction could be exploited 
further to test the possibility that the series can take simpler forms, as compared to a 
brute force expansions to inhnite order. We leave this interesting question for future 
investigation. 


f.3. Currents found from the Rayleigh Sehrodinger (locator) expansion 

A natural question that arises is the relationship between the currents found above and 
those found from a brute force expansion of the projection operators of the Anderson 
model in powers of the coupling constant y. The model has a formal single particle 
eigenfunction expansion in the form 

1^(2/)) = X^“ofc(2/)4lO)i (32) 

k 

with an initial condition localized say at the site 0 as t(ofc(0) = Sko- The projector 
Q = |\l/(|/))(\l/(|/)| can be expanded in a series in y 

Q = Yi “Oj“otc]ct = P(0) + yP(l) + y‘‘P(2) + ... 


(33) 
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so that the basic expansion of the wave functions in a Rayleigh Schrodinger (RS) series 
in y generates the conserved currents. We can use the standard result in text books 0 
to write a perturbative expansion for the state at site 0 with standard normalization to 


'Uqq — 1 clS 

l^(l/)) = 4|0>+ X^wofc4|0>> (34) 


with a power series expansion for uofc 


tok 

uok = -y - 

€o ~ 


1^0 


_ tkltio _ 

(eo - efc)(eo - e/) 


2 fooffcO 
(eo - efc)2 


Oiy^). (35) 


Using this expansion, we may generate the series Eq. (jSHl), the result is explicitly stated 
below in Eq. flHSll . Erom this series we can verify to second order, that this series differs 
from that in the standard gauge Eq. (TTnil by specific gauge terms. The advantage of 
Eq. ffT^ is that this gauge invariance is manifest in the construction by the nested 
commutators. On the other hand, Eq. (13311351) . corresponds to a particular gauge picked 
out by the R-S method, and the currents found here are some linear combinations of 
the ones in Eq. (fT9l) as in Eq. fl23l) . 

It seems to us that the series in Eq. (in possesses an essential simplicity relative 
to the Rayleigh-Schrodinger series Eqs. (13311351) . The R-S perturbation expansion 
simultaneously determines the energy eigenvalue, and for this purpose very specific 
gauge terms are needed. On the other hand all terms in Eq. (TTni) are generated by 
completely off diagonal terms, those terms that avoid multiple visits to any site. This 
leads to simpler recursion relations, as in Eq. (1^ . relative to the RS series. Eor this 
reason our numerical work in this paper uses the series in Eq. (fT9l) . 


4-4- Locator expansion for Type-1 Hamiltonians 

The Rayleigh-Schrodinger series can also be constructed for Type-1 Hamiltonians using 
the exact eigenstates IE) with eigenvalues E as given in Eq. ([71) and Eq. (JHl). The 
projector \E){E\ = (gZ'Jxg-e-) expanded in y as shown in Eq. In the 

limit r/ —7- 0, the roots of Eq. (jHl) tend to e*. We take the root E ^ cq to obtain Qq 
the conserved charge corresponding to site 0 calculated using the Rayleigh-Schrodinger 
gauge. Other roots yield other Qj. Expanding Eq. (jH]) for E in y near E = cq, we get 

^ = €o - |/7o + l/^7o X] +0{y^). (36) 

~ 


^ For e.g. see Eq. (5.1.44) Modern Quantum Mechanics^ J J Sakurai (Pearson Education 1994)) 
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Since, the projector diverges in y ^ 0 limit, we define our conserved charge as 
Qo = \E){E\ to make it well behaved. Qq is given by, 


~ -E" — €o 7i(coCj + cjco) 

<3o = no + ^— 2^ 


+ 


.VO 

{E - 


E-Cj 


7o 


E 

*jV0 


lajclcj 


{E-e,)iE-e,) 


(37) 


Then, replacing - > ?/^7o ^-nd then E ^ cq, we have obtained Qo as a combination 

of Qi (see Eq. IS]) as follows, 


Qq — Qci — y 


E 

z/0 


liQi + liQci 


Co 


+ 0(?/3 


(38) 


Other Qfc can be obtained with the replacement 0 —?• fc. Unlike Qfc, there is no indication 
of the series truncating at any finite order for Q].. 


4 . 5 . Convergence of the power series 

The conserved charges constructed above depend on the microscopic parameters of the 
Hamiltonian, i.e. the hopping and on-site energies. As we shall show later, the same 
Hamiltonian can have a localized and delocalized phase depending on the values of 
these parameters. It is thus important to understand if and how the conserved charges 
themselves differ in the two phases. More precisely, how do the conservation laws 
“know” whether a particular choice of microscopic parameters produces a localized or 
delocalized phase? 

The answer has to do with their convergence since they are expressed as power 
series in the microscopic parameters and particle operators. We thus need to state in 
what sense the power series are convergent. A reasonable condition for convergence 
is a sufficiently rapid decay of the coefficients rjff with increasing m. However, this is 
complicated by the fact that there are energy difference denominators in the coefficients 
rfff that can cause them to blow up when the on-site energies at two different sites are 
equal. To avoid this, we restrict ourselves to a particular type of disorder that may 
be termed ” non-resonant”. By this we mean any ensemble of e,, which shows “level 
repulsion”, i.e. the probability of finding e* very close to each other is very small. 

From the random matrix theory, we know that the eigenvalues of a generic matrix 
display level repulsion in their eigenvalues of various degree, the Gaussian Orthogonal 
Ensemble (GOE) [SO] of real symmetric matrices has the least level repulsion. This 
condition ensures that perturbative resonances from small denominators, that would 
otherwise cause individual terms in the expansions of the conserved charges to diverge, 
are prohibited. This choice is similar to the one involving limited level attraction recently 
adopted in the context of many-body localization Ell. 

We have verified that this distribution of onsite energies gives us localization (as 
indicated from a calculation of the participation ratio) immediately upon switching on 
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Figure 3. Plot indicating convergence of conserved charges [see Eqs. m and 
of the Anderson model ([T 8 |) for N = 500. \r]\ represents a typical m-th coefficient 77 ^ 
averaged over a distribution of the on-site disorder e^, see the end of Sect. 14.51 The 
plot shows the logarithm of the average as a function of m. are drawn from the 
eigenvalues of real symmetric matrices whose elements are Gaussian random variables 
of variance a =0.1, 0.25 and 0.4, and we set yt = 1. 


the hopping term. Thus, this particular choice of onsite energies, which is of great 
convenience from the point of view of calculations, is also not nnphysical. The on-site 
energies are drawn from the eigenvalues of a real symmetric matrices whose elements 
are taken from a Gaussian random distribution with fixed variance. The eigenvalues of 
these matrices are assigned randomly to different sites. Different random assignments 
then constitute different realizations of disorder, which can then be averaged over to 
check for convergence. The result of this procedure is shown in Fig. [3l where are 
drawn from the eigenvalues of real symmetric matrices whose elements are taken from 
a Gaussian distribution of variance ( 7 = 0 . 1 , 0.25 and 0.4. It can be seen that the 77 ^ 
decrease rapidly with increasing order of power series m indicating convergence. We 
have also checked the convergence of the power series for drawn from the eigenvalues 
of non-integrable t — t' — V model, which also follow a GOE distribution [33l [5^ . 

Since 77^ contain more than one term for each 772, we checked the convergence of a 
typical term, which is of the form 7 -w- —7 - 7 . Recall that the 772 ^^ order 

^bi)\^a2 ^£>2 / • • • • Vctm ^bm) 

term in the calculation of Qo involves sites with labels between N — {m — 1) and m as 
can be seen from Fig. [ 2 l Thus, the only values of e* involved are are those chosen from 
[ejv-(m-i)) Cm] (co Is at the center) such that 7 ^ e;,., Wi and max |aj — 6 j| = m 

As the aim of this work is to construct conserved charges in localized systems, it 
is legitimate to ask whether this slightly non-standard choice of disorder distribution 
produces localization. We have verified this through numerical exact diagonalization by 
calculating the PR. We find that the PR for different eigenstates is indeed close to zero 
for systems of size N = 500 as shown in Fig. 01 consistent with localization. We thus 
conclude that our model with on-site energies taken from a GOE distribution does indeed 
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Figure 4. PR of eigenstates of the Anderson model (p!8|) for N = 500 numbered in 
ascending order according to the energy levels. On-site disorder is drawn from the 
eigenvalues of real symmetric matrices whose elements are Gaussian random variables 
of variance a =0.1, 0.25 and 0.4, and we set yt = 1. Blue dashed line corresponds to 
the typical value of PR in delocalized phase. 



Figure 5. Conserved charges for Aubrey-Andre model converge for h > 1 (localized 
phase) and diverge for h < 1 (delocalized phase). Here \r]\ represents a typical m-th 
coefficient in Eq. (|20|) [see the end of Sect. 14. 5j . N = 900 and (3 = The plot 

shows log 17^ I as a function of m. 


produce a localized phase. A similar exercise to construct the conservation laws for the 
above model has been carried out in Ref. [53]. In that work too, the conserved charges 
have been constructed as infinite operator series but whose coefficients correspond to 
the amplitudes of a particle to be on the sites of a square lattice whose sides are the 
physical one dimensional lattice. The recursion relation obtained is between conserved 
charges on different sites and the convergence of the series is assumed to follow from 
the exponential decay of the eigenfunctions of the Hamiltonian. In our calculations, we 
construct the conserved charges directly in terms of the microscopic parameters of the 
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Hamiltonian and our convergence criterion is not based on any assumption about the 
nature of the eigenstates of the Hamiltonian. In fact, as we show in the next section, the 
convergence of the series for the conserved charges can be used to identify the delocalized 
and localized phases instead of the eigenfunctions. 

5. Aubry-Andre model 

Having constructed the conserved charges for a model with finite-range hopping and 
defined a condition for convergence of the power series for them, we can further 
investigate the meaning of our convergence criterion. In particular, since our goal is 
to identify the validity of our construction of the conservation laws with the presence 
of localization, the power series should fail to converge according to our criterion in a 
delocalized phase. 

We thns require a non-interacting model with disorder in one dimension which 
has a delocalized phase. While any model with finite-range hopping and an on-site 
random potential in one dimension always produces localization BE], a quasi-periodic 
potential can produce localized and delocalized phases. Such a model is the Aubry- 
Andre model given by the Hamiltonian 

H cos(2/37rj)c]9 - ^ ^^(cjcj+i + h.c.), (39) 

j j 

where (3 is an irrational number. The parameter h can be tnned to effect a transition 
from a localized phase (for h > 1) to a delocalized phase (for h < 1) [12]. We note that 
this model is usually studied with an additional term that introduces a p-wave pairing 
gap |5lj, but we set it equal to zero for our analysis. 

The localized phase here is one in which all single particle states are localized and 
similarly all single particle states are delocalized in the delocalized phase. The transition 
between these phases happens at h = 1. Since the Hamiltonian in Eq. fl3^ is also of 
the form inai, we can use the expressions obtained for the riff in the previous section 
to construct the conserved charges. These will now depend on the parameter h (i.e. 
y —7- (2h)“^ in the previous section) and if the criterion for convergence postulated by 
us is a valid one to detect localization, we shonld observe the power series to converge 
in the localized phase [h > 1) and diverge in the delocalized phase [h < 1). This is 
indeed the case as we see e.g. from Fig. [5l which shows that a typical matrix element 
of goes to zero quite rapidly with increasing m for /i > 1 but diverges for /i < 1. 
Thus, we have established that our convergence criterion is valid for identifying the 
localization-delocalization transition. 

6. Interactions 

We now turn to systems with interactions. The simplest way to introduce interactions 
to models we studied here is through a nearest neighbor density-density term. Let us. 
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Integrals of motion for one-dimensional Anderson localized systems 
for example, add such a term to Eq. (ITHD . 

H CiUi - ty^{c\ci+i + h.c.) + V^ n^ni+i 

i i i 

= Hq + VSH, (40) 

where we redefined Hq as compared to Eq. ing. 

We assume that the particles here are spineless fermions. It is tempting to try a 
construction of the conserved charges starting from a zeroth order Hamiltonian that 
combines the on-site and interaction terms since they commute with each other and 
their eigenstates are localized at every site. However, the interaction term is quartic 
in creation and annihilation operators and so the conserved charges can no longer be 
assumed to be power series in the hopping with each term quadratic in the creation and 
annihilation operators. Such an assumption leads to no solution for the coefficients since 
the commutators keep producing terms with increasingly longer trails of creation and 
annihilation operators as one goes to higher orders in the hopping. A more profitable 
exercise is to try to obtain the conserved charges as power series in the hopping but 
only to the first order in the interaction. While these are not exact, they offer a 
reasonable approximation in the limit of small interaction strength. Weak interactions 
typically should not destroy the localization present in the non-interacting limit and 
thus conserved charges should continue to exist. 

We know from our previous calculation that the operator of the form Qq = 
no-\-'f2ijm commutes with Hq. Let us now define a new operator Q = Qo-\-VSQ 

to linear order in V and calculate the commutator. 

[Q, H] = [Qq + V5Q, Hq + VSH] 

= V {[SQ, Hq] + [Qq, SH]) + 0{V^). (41) 

We choose SQ such that [SQ, Hq] -|- [Qo, SH] = 0, so that Q and H commute to 0{V). 
We assume the form SQ = Note this is quartic in the creation and 

annihilation operators since the interaction term is as well. Thus, 

[SQ,Hq]^[Qq,SH] = {) 
^y'^{.'n^c\cknk+i - rikAvnk+i + rif^k+\'^kc\ck^i 

kim 

Vk-[-l,i^k^k-^l^i) ^ "^rst—lv "^rstv+l 

rstv 

^ ^ "^rstvi^^v ^J• 

rstv 

We now assume that 'iprstv can be written as a power series in y, i.e. 'iprstv = 
' Equating the coefficients at different orders of one can in principle 
obtain in terms of the for the case with E = 0. fn fact, it can be seen that at 
a given order m, the are linear combinations of the and the A^^Ai> One can 
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also impose constraints arising from the anti-commutation of the fermionic operators, 
the Hermitian nature of the conservation laws and the number of non-zero components 
of the to severely constrain the number of non-zero components of 

Let us, for example, derive 5Q to the first order in r/, i.e. we set m = 1. We have 

'^{vlAckrik+i - rilAkCm+i + - r}l^i^^kcl+iCi) 

ki 

= -\- tt — tg + er)clcsclcy. 

rstv 


Since, only rjl^j^i and are non-zero, the non-zero A^^"^ are given by the following 
equations: 

^ k-\-i-,k^k-\-i-,k-\-^ ^ k,k-\-^,k-\-i-,k-\-^ ^k-\- 1,fc 

-;; r 

€fc+l — 

^k,k—l,k+l,k+l _ ^k—l,k,k+l,k+l _ 'Ik—l,k 


Cfc-l — Cfc 
1 


j^k,k,k+l,k _ j^k,k,k,k+l _ ^k,k+l 

Cfc — Cfc+l 

j^k,k,k+2,k+l _ j^k,k,k+l,k+2 _ ^fc+2,fc+l 

^k+2 ~ ^k+1 


The corresponding expression for 5Q to order y is 
5Q = yy^ AY^^'clcgclcy = yV^ 

rstv k 

t \ ^k —1 k / i t \ 

+CfcCfc+i)nfc+i H-^—(c^_iCfc -|- c^Cfc_i)nfc+i 


(c^ c, 
^ ^ \^k+l^k 

l^k +1 — ^k 


^k -1 ~ ^k 

+ UkicAck + 4cfc+i) + 


^k ~ 


^k +2 ~ 

(4+2^^+! + cl+lCk+ 2 ) 


Other approaches to construct conservation laws for interacting systems have been 
proposed including a recent one where the interacting problem is mapped onto a non- 
Hermitian problem on a lattice in operator space [S2]- A convergence criterion for 
the resultant series based on the operator norm is then used to identify localized and 
delocalized phases. 


7. Conclusions and discussion 

Inspired by the Type-1 Hamiltonian system, we have demonstrated a scheme to obtain 
the conserved charges for non-interacting disordered models displaying localization in 
one dimension. One of our motivation was an observation of similarities between 
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localized and integrable systems, such as the absence of level repulsion and the absence 
of thermalization. Our conserved charges are exhibited as a power series in the hopping, 
and using a suitable convergence criterion, we show that the convergence (or divergence) 
of conserved charges tracks the presence (or absence) of localization. An interesting issue 
of “gauge dependence” of the conserved charges is unearthed and explored. It is shown 
that a full understanding of the gauge dependence leads to considerable simplifications 
of the charges in some cases. On the other hand, straightforward Rayleigh Schrodinger 
perturbation theory or equivalent schemes, commit one to a particular gauge that is 
often inconvenient. 

This work provides a novel link between the concepts of localization and 
integrability. Our results hold within the context of the 1-d Anderson model, where all 
states are localized, and the Andre-Aubry model, where (all) states undergo a transition 
tuned by a coupling constant. It is not immediately obvious how to extend these 
results to a higher dimensional Anderson model with a mobility edge separating the 
two classes of states. The Aubry-Andre model exhibits an interesting kind of duality 
which allows the localized and delocalized phases to be mapped onto each other with 
the roles of the hopping and onsite potential exchanged. The duality transformation 
is expressed in terms of new fermonic operators given by cj. = exp(z27rfc/3)c„, 

which are eigenstates of the momentum operator with eigenvalue: k = kF^-i mod F„, 
where Fn is the n-th Fibonacci number and L = Fn [SSI EE] • In terms of these fermionic 
operators the Hamiltonian fl3^ becomes 

^ cos{2P7rk)n-k - - ^(4cfc+i + h.c.). (42) 

k k 

The Hamiltonian satisfies the duality relation: H{h)/h = H{l/h). We have shown 
that for the Aubrey-Andre Hamiltonian written in real space, one can construct set of 
conserved charges that converge for 0 < h < 1. Because of the duality of the model one 
can construct similar conserved charges in terms of cj. and 4 . The power series of these 
charges converge when 0 < l/h < 1 and both sets of charges diverge at /i = 1. Thus, 
the duality of the model allows us to explicitly construct conservation charges in one 
phase given that they exist in the other. 

This can be better understood by noting that localization is a basis dependent 
concept. We have been using localization (as is the standard practice) to mean 
localization in real space. To obtain the conserved charges for such a localized phase, 
we start from a Hamiltonian whose eigenstates are perfectly localized in real space 
and then add terms perturbatively in the hopping. Similarly, the delocalized phase 
of the Aubry-Andre model is localized in momentum space and one can then obtain 
its conserved charges by starting with a Hamiltonian perfectly localized in momentum 
space (tight binding model) and then add terms perturbatively in the on-site potential. 
This is the essence of the duality outlined above. Thus, the conserved charges also 
carry labels indicating the space (real or momentum) where the system is localized. 
What is important though is that once the basis in which the system is localized is 
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identified and the conserved charges are constructed accordingly, they are sensitive to the 
onset of delocalization in that basis and can be used to locate localization-delocalization 
transitions. 

The importance of the basis can be further understood when one compares 
the behavior hard-core bosons with that of spinless fermions in the Aubry-Andre 
model [SZlEH!. The duality between the localized and delocalized phases is destroyed 
for hard-core bosons. As a result, the relaxation of real space local observables in the 
localized phase is different from their conjugates in momentum space in the delocalized 
phase. This feature is absent for spinless fermions where the duality holds and as a 
consequence, conserved charges of the type derived in this work exist in both phases. 

While it is only possible to construct these charges to lowest order in the interaction 
using our procedure, their fate upon the introduction of interactions can in principle be 
investigated numerically, which we defer to a future workEI 
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